Non-adiabaticity of a many-particle Landau-Zener problem: forward and backward 

sweeps 
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We consider nonlinear dynamics of an interacting many-particle system with a slowly changing 
parameter: a many-body generalization of the Landau-Zener model. Adiabaticity is destroyed when 
the parameter crosses a critical value. Applying semiclassical analysis based on concepts of classical 
adiabatic invariants, we derive a formula which accurately describes particle distributions in the 
Hilbert space at wide range of parameters and initial conditions of the system. We found striking 
universal features in the particle distributions which can be probed in an experiment on Feshbach 
resonance passage or a cavity QED experiment. The dynamics is found to be crucially dependent 
on the direction of the sweep. 



Destruction of adiabaticity is a central issue in both 
quantum and classical mechanics. For a single-particle 
quantum system, very often nonadiabatic dynamics can 
be described within the exactly solvable Landau-Zener 
model (LZM) Q, where a probability of transition from 
an initially occupied instantaneous ground state level to 
the excited one is exponentially small in the sweeping 
rate parameter. Experiments with ultracold quantum 
gases may involve macroscopically large numbers of 
particles, making semiclassical (SO) treatments justified 

iiiiail. 

Nonadiabatic dynamics in SC models of many-particle 
systems has been intensively discussed recently 0, [(| 0, 
M,M- I n many treatments, it was found that exponential 
LZM-typc behaviour for the transition probabilities is re- 
placed with power-laws of the sweeping rates. However, 
a universal and accurate method for describing nonadia- 
batic SC dynamics of such many-particle systems is still 
lacking. Here we present a method based on modelling of 
a quantum system by an ensemble of classical trajectories 
and accurate description of its dynamics in a vicinity of a 
bifurcation, which predicts universal features of the par- 
ticles distributions amenable for an experimental probe. 



We consider a many-particle LZM: time-dependent 
Dicke model, which has numerous applications in quan- 
tum and matter- wave optics @, H, H, Mill- The model 




FIG. 1: The time evolution of the distribution of the rescaled 
boson number n at a "forward" sweep through the resonance, 
(a) quantum calculation with N=500. (b) Corresponding en- 
semble of classical trajctories; sweeping rate is e — 0.3 
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ftb + l^-S' + ^={PS- + bS+), (1) 



are 



where is the coupling strength, = S x ± iS y 

spin operators, b' and b are creation and destruction op- 
erators of a bosonic mode, j(t) = ±2et is detuning, and 
e is the sweeping rate of the bosonic mode energy. The 
spin value S is macroscopically large S = N/2 ^> 1; usu- 
ally the physical origin of the effective spin variable S 
is a collection of two-level systems (spin-1/2 particles). 
Among many possible applications, let us mention the 
Feshbach resonance passage [l3j], dynamics of molecular 
nanomagnets [TJ], and cavity QED with Bose- Einstein 
condensates (BEC) HGJ]. 

With N = 1, one recovers the standard Landau- 
Zener model. In the context of a Feshbach reso- 
nance passage in a Fermi gas, an equivalent realization 
of the Hamiltonian (JTJ is -ff = d> J2iLi (™*T + + 
X^fcLi f^c,'J.CiT + h.c?j , where c; iCr and b are the 
fermion and boson annihilation operators, respectively, 
hia = cJ^Cjo-, and a =|, J.. Coupled atom-molecular BEC 
are described by a similar model, whose dynamics in the 
limit of large N becomes equivalent to ([1]) with the re- 
placement 7 — ► —7 [il [To|. I.e., in the thermodynamic 
limit of the degenerate model ([T]) association of Fermi 
atoms becomes equivalent to dissociation of molecular 
BEC, and vice versa. We therefore discuss two driving 
scenarios here: "forward" and " backward" sweeps. In the 
"forward" sweep, starting in the distant past with some 
small initial number of bosons Nb{t)\_ = (&'' 6) (t) | = 
N- = n N, we want to calculate the final number of 
bosons Nb(t)\+oo = nfiV and its distribution P(Nb) as 
a function of the sweeping rate e and the initial bosonic 
fraction n_. To be more specific, we start in a number 
state of ©: |V_oo) = \N-)\S,S,) = \N_) |f , f - JV_), 
where N_ can be zero in case we start at the ground state. 
We consider the sector S — N/2 for clarity, but other 
values of S can be treated analogously (Hilbert space of 
([1]) is decoupled on sectors with definite total spin, or 
"cooperation number", S). Matrix elements of (fl]) have 
= -7(t)5„,„' + n8 n , n , +1 ^N - n'/VN + 



the form H n 
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n'S n ^ n i_i\/N — n/y/N (we set g = 1 for convenience in 
this Letter), therefore the model can be referred to the 
class of generalized Landau-Zener models [l6[ in the case 
of linear driving ^(t) ~ t. In the second dynamical sce- 
nario ("backward"), we start in the ground state of (fT]) 
at large positive value of 7 and make an " inverse" sweep 
to large negative 7 (see also our supplementary Brief Re- 
port [lOj]). This scenario is relevant to association of Bosc 
atoms in molecular BEC, studied e.g. in [6|. 

The model and its equivalent realizations have been 
studied already using various techniques. Diagrammatic 
methods employed in 0] work well for small N but do 
not allow to get close to the SC limit for large N. An- 
other approach is a SC treatment based on classical adi- 
abatic invariants @, B 0|- At \t\ = 00 and N — » 00 
the relative number of bosons N^/N corresponds to a 
classical action I of an effective Hamiltonian system (see 
e.g. Eq.([2])), so that for large but finite N the problem 
can be mapped to the calculation of a change of classi- 
cal actions of a properly prepared ensemble of trajecto- 
ries. It is inspiring that calculations with ensembles of 
classical trajectories (defined below) reproduce full quan- 
tum calculations very well (see Fig[T]) . The classical limit 
of ([1]) can be obtained by introducing Schwinger bosons 
d,/3 : S + — a' (3, S~ = (3* a, S z = (n a — ng)/2, 
and using the c — number formalism 0, Q. The re- 
sulting classical system has 3 pairs of variables (t1q!,/3,6 
and conjugated phases) and 2 integrals of motion (Af a = 
n a + np, L = np — rib). We reduce the system to 1- 
degree of freedom Hamiltonian: 



H = --yn - 2y / n(n + L)(N a - L - n)cos0, (2) 

where n € [0, 1] corresponds to the rescaled number of 
bosons Nb/N, and <fi is the canonically conjugated phase. 
We employ the Truncated Wigner Approximation [8|, [TtJ 
and prepare our initial ensemble of classical phase points 
using Wigner transform of the quantum ground state at 
large negative 7. Far from the resonance, the quan- 
tum system is effectively decoupled on 3 oscillators, 2 
of them being in the ground state, and the third one in 
the high- lying Fock state \N). With a good accuracy, 
W(nb !a ,p) ~ exp(— 2nb — 2np)5(n a — 1). It means ini- 
tial phase <f> is uniformly distributed on (0, 2tt), and be- 
sides exponential distribution of initial n we have to take 
into account distribution of the parameter L. Consider 
firstly the case L — 0, then the Hamiltonian reduces to 
H = —771 — 2n»Vl — ncos(j>. Note that by means of a 
trivial change of variable sJ18l | , the Hamiltonian becomes 
the same as analyzed in [9j, and can be derived from SC 
Bloch equations of The classical phase space of the 
Hamiltonian ([2]) at fixed values of the parameter 7 is de- 
scribed in (9I llOj. If 7 < —2, there is only one stable 
elliptic point on the phase portrait (Fig. 2a). At 7 = — 2, 
a bifurcation takes place. There are two saddle points 
at n = 0,cos</> = —7/2, and a newborn elliptic point at 
= 0. The trajectory connecting these two saddles (the 
separatrix) separates rotations and oscillating motions. 
At large positive 7, again there is only one elliptic sta- 
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FIG. 2: Phase portraits of ©. From (a) to (c): 7 = 
— 10,-1.4,20 correspondingly. Saddle points are denoted by 
asterisks, the bold line is the separatrix. Shaded areas illus- 
trate definitions of the classical actions (see text). 



tionary point at <j> = 0, and n close to 1. The classical ac- 
tion is defined as in [8|, [9| and is shown graphically in Fig. 
2: shaded areas, divided by 27r. At 7 = —00, the action 
coincides with n: I = n, and canonical "angle" variable 
(p coincides with <fi. At 7 = +00, we have I = 1 — n. We 
are interested in small initial actions, corresponding to 
the ground state of the quantum system. 

We split initial classical ensemble on slices with equal 
actions L_ (let us call such an ensemble Aj_), and an- 
alyze dynamics of each slice in detail. Initial actions 
i_ of a classical ensemble in the case of a fully polar- 
ized (n_ = 0) initial quantum state are of order due 
to quantum fluctuations. In case there is some non- 
zero initial population of the bosonic mode, 7_ can be 
much larger than jj. There are three small parame- 
ters in the model: the sweeping rate e, and the ini- 
tial action i_. Naturally, 4 < e, Still, freedom of 
choosing a value of the ratio J_/e remains. Assuming 
i_ 3> e and e, l/N, L<1, which corresponds to consid- 
erable initial population of the bosonic mode, Ref.0 ap- 
plied an approach based on the separatrix crossing theory 
[l9t I20I I23I ]. From the experimental point of view, how- 
ever, it is important to analyze also cases of very small 
initial actions. 

Our method for the regime of small initial actions is 
as follows. To calculate the deviation from adiabatic- 
ity, we note that most of the change of the classical 
action of a phase point happens as it travels near the 
separatrix and, especially, near the saddle points that 
arise during the bifurcation as 7 reaches — 2_ § (in our 
new variables, the arising saddle point is located in the 
origin). Close to the saddle point n is small and one 
can expand y/\ — n in series. Let us introduce variables 
P = Y = 2y / nsin^. Near the bifurcation, 

at 7 ps —2, we therefore get an effective Hamiltonian 
H = -^(7/2 + 1) - £(7/2 - 1) + fj, where higher 
order terms in Y and P have been neglected. We neglect 
also the time-dependence of the coefficient 7/2 — 1 of the 
second term. Then, shifting origin of time, we obtain the 



Hamiltonian 



-et 



Y 2 + P 4 /16. After simple rescal- 
= 2 2 / 3 F',t = t'/2 1 /3 j H = 2 7 / 3 H\ 



ings P = 2 4 / J P', Y 
the Hamiltonian becomes H — — -^-ei+^- + -^- (primes 
over new variables omitted). We replace also momenta 
and coordinates: P — ► Y. Y — ► —P. Let us now in- 
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FIG. 3: (a) Change in the action AI of classical trajectories 
from the ensemble Ai_ as a function of the initial phase. 
The change AI describes the deviation from adiabaticity. We 
shift all curves horizontally so their maxima are at 7r/2. We 
found the change in the action is highly phase-dependent, 
(b) Average change of the adiabatic invariant {AI)^! as a 
function of the sweep rate e. Lines are the theoretical curves 
(Eq[SJ). Squares, triangles are numerics for J_ = icr 3 ,i(n 4 , 
respectively, (c) Number of created bosons from quantum 
calculations with N = 250 and N = 500 (dots); from the 
Formula [8] obtained by averaging over the initial distribution 
W ~ exp[— 2nb — Znp] (solid lines), (d) Dispersion of number 
of bosons for N=250; dots: quantum calculations, solid line: 
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the formula ©, i.e. N 2 (n 2 - n' 2 



troduce a rescaling transformation that makes the es- 
sential mathematics of the problem as clear as possible: 
Y = e^ 3 Y, P = e 2 / 3 P, t = e" 1 /^, 77 = e 4 / 3 H. 
The Hamiltonian becomes (omitting tildes over new vari- 
ables) H — — s^- + -^- + ^-. This Hamiltonian does not 
have a small parameter any more, and the loss of adia- 
baticity is evident. An important property of the bifurca- 
tion we are considering is that the effective Hamiltonian 
leads to the second Painleve equation (PII) 



<PY 
ds 2 



= sY- 2Y i 



Asymptotics of PII were investigated by Its and Kapaev 
|22| (see also |24l) using a method of isomonodromic de- 
format ions|25[ . At s — * — oo the asymptotic solution to 

© is m mr 



Y(s) = „(-.s) -sin ( |(-s) 3 / 2 + ^a 2 ln(-s) + . 



and in the limit s — > +oo it is 



VM ^/£±p(2 S )-icos(^ S 3 / 2 -^ 2 ln( S ) 



where (a,4>) and (p,9) are the "action-angle" variables 
characterizing the solutions in the limits s — > ±oo. As 



s — > ±oo, the adiabatic invariant I p of equation ([3]) ap- 
proaches the quantities 7~ or 7+ which are defined (to 

the lowest order terms) as I~ — 7+ = The 

jump in the adiabatic invariant AI P = 27+ — I~ can 
be found from general relations between p 2 and a 2 as 
p = Ve 

[221 . |24|). Returning back to the original variables and 
the Hamiltonian ©, we get the Formula: 



J P 2tt 2|lm(p)|' 



l e U(I P )-i<t>] ( se e 




- 1 ln(2sinO 

(4) 

where for the ensemble Ai_ 7r£ = (/(7_) + the func- 
tion /(7_) is not important for our discussion (see [13]); £ 
is a quasi-random variable uniformly distributed on (0,1). 
The formula predicts the average change in the action to 
be 



7_ 



In i / exp 



tt7_ 



(5) 



which means that the final number of bosons is n/ = 
1 — 7_ — (A7)_4 7 . It predicts also the final distribution 

P(A7) = P(l - 7_ - n s ). All moments A7^ f = ((A7 - 
(A7)) fc )_4 7 are easy to calculate, for instance M\ = 

(f) 2 / 1 ^ln 2 (2sin<) = 4. 

Note that the phase-dependent part in (j4|) exactly co- 
incides with the result of Q , obtained under the different 
conditions e « L < 1. Moreover, if 7_ ^> e, we recover 
this result of Q, which in the present variables is 



A7 



2c 



ln(2sin 7r£), 



(6) 



where the pseudophase £ € (0,1) is a quasi-random vari- 
able. Such a change in the action has zero mean value, 
nevertheless it introduces spreading in particle distribu- 
tion since (A7 2 ) ~ e 2 . When 7_ <C e, we have a quali- 
tatively different result which resembles that of H (i.e. 



(3) 7+ = eln7_/7r, with 7_ = i) 



A7= -e(-ln( — ) + - ln(2sin7r£) 
\7T V e / 7T 



(7) 



Formula ((4]) can be used over a wide range of val- 
ues of e and 7_ (we require 7 « 1 and e|ln7_| < 1 
). Qualitatively, it is important that the final distribu- 
tions are determined not only by the average change in 
the action ((5|), but also by the phase-dependent part ([6]), 
which is therefore important in amplification of quantum 
fluctuations. This profile has striking universality: we 
found that several other infinite-range models, e.g. Lip- 
kin model, two-mode BEC with attractive interaction, 
have analogous behaviour during linear sweep of a con- 
trol parameter through a critical value. The comparison 
between the numerical and analytical results for 7 = 



4 



are given in Figs. [3^,b, where predictions of the Eq.Q 
start to deviate from the classical numerics only at large 
e (such that — ^ In ^f- ~ 1). Then, we take into account 
distribution of L. We consider a slice of total distribution 
with equal values of x = n + 4 and uniform distribution 
of L G (— 2x, 2x); averaging over <fi and L within the slice, 
we found that the formula([5]) (derived for the "central" 
point of the slice, i.e. L = 0,x = J_) acquire addi- 
tional coefficient 2 inside the logarithm: I_ — > 2i_. The 
amplitude of the phase-dependent part is also modified, 
although it retains its characteristic form shown in Fig. 
3a; the final result of averaging over the 3-dimensional 
distribution (in <j>, 1_ , L) of the initial ensemble of phase 
points is: 

n=l-~ In— + 7 , n 2 -n 2 = -, 8 
n \ n J 6 

where 7 is the Euler constant. Comparison with quantum 
calculations are given in Figs. 3c, d. 

Let us now briefly discuss the inverse sweep (see also 
[Iol|). whose properties can also be derived from PII equa- 
tion. The main feature of the inverse sweep is phase- 
independence of the change in the action in the limit of 
small I—. Deviation from adiabaticity at slow sweeps 
strictly follows a linear power-law. The coefficient was 
estimated in [6j to be equal to ^ w 0.21. The asymp- 
totically exact value of the coefficient was found from 
mapping to PII in 10]: ^ w 0.2206. The most appeal- 



ing physical system to implement this process is photo- 
or magneto- association of Bose condensates, which was 
shown to be well described by a two-mode mean-field 
model. "Forward" sweep would correspond to dissocia- 
tion of BEC. Experimentally, studying distributions after 
a sweep could be achieved by making many runs at each 
particular sweeping rate. 

In summary, we have found novel universal features in 
dynamics of a many-particle LZM using a method that 
complements the separatrix crossing theory applied to 
that model in [9[. The results are highly relevant for ac- 
curately describing the Fcshbach resonance passage in ul- 
tracold Fermi and Bose gases and may motivate further 
experimental activity on that theme. Furthermore, we 
believe our method will have important applications in 
the field of dynamics of quantum phase transitions (QPT) 
[26Tj . When an underlying quantum system experiences a 
second order QPT, its mean-field counterpart often ex- 
periences a Hamiltonian pitchfork bifurcation. Passage 
through such a bifurcation leads to PII. 
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